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Abstract 


An electronic throttle is a DC servo drive which positions the throttle plate, thus providing drive-by-wire control of engine torque. 
This paper presents an electronic throttle control strategy consisting of a PID controller, and nonlinear friction and limp-home 
compensators. The emphasis is on the development of an adaptive control strategy, which is aimed to enhance the control strategy 
robustness with respect to process parameter variations, caused by production deviations, variations of external conditions, and 
aging. The adaptive strategy consists of auto-tuning and self-tuning algorithms. The auto-tuner provides automatic tuning of the 
control strategy parameters during vehicle assembly, or uses specific drive modes during normal use (e.g. each time the engine is 
turned off). The self-tuning algorithms are based on the permanent, on-line estimation of those process parameters which can vary 
within a single engine run (e.g. the DC motor armature resistance). The presented adaptive control strategy has been verified 


experimentally. 
© 2005 Elsevier Ltd. All rights reserved. 
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1. Introduction 


In conventional vehicles the driver gas pedal is 
mechanically linked to the throttle plate. This system 
is being replaced in modern vehicles by an electronic 
throttle, where the link between the gas pedal and the 
throttle plate is realized by means of a DC servomotor. 
In this way, the engine control unit can correct the 
throttle position reference value for specific engine 
operating modes, thus improving drivability, fuel 
economy, and emissions, and also providing the 
implementation of engine-based vehicle dynamics con- 
trol systems including traction control (Huber, Lieben- 
roth-Leden, Maisch, & Reppich, 1991). The electronic 
throttle consists of an electronic throttle body, a 
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chopper, and a position control strategy. A photograph 
of an electronic throttle body is shown in Fig. la, and its 
internal structure is depicted in the lower right part of 
Fig. 1b. The electronic throttle servo system does not 
usually include an inner current controller. The throttle 
motion is constrained by a dual return spring which 
returns the throttle in the so-called limp-home (LH) 
position in the case of power supply failure. 

It has been demonstrated in (Deur, Pavković, Perić, 
Jansz, & Hrovat, 2004) that DC drive transmission 
friction and the dual return spring nonlinearity at the 
LH position significantly affect the performance of the 
electronic throttle control system. A nonlinear control 
strategy is proposed in the same reference in order to 
compensate for friction and LH nonlinear effects, thus 
providing a linear system-like behavior for a wide range 
of throttle operation. The step response of control 
system is characterized by the settling time of approxi- 
mately 70 ms and the steady-state accuracy better than 
0.1°. The control strategy has been tuned based on the 
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Fig. 1. Photograph of electronic throttle body (a), and principal schematic of adaptive electronic throttle control system (b). 


results of off-line experimental identification of electro- 
nic throttle process model (see Pavkovic, Deur, Jansz, & 
Perić, 2003). 

It can be expected that the performance of an 
electronic throttle with constant controller parameters 
will deteriorate in the presence of process parameter 
variations. There are three sources of process parameter 
variations (see Pavkovic & Deur, 2002; Hashimoto, 
Ishiguro, Yasui, & Akazaki, 2003): 


(a) Production deviations. Different ETBs from the same 
production series can have different parameters. 
Static curves u(@) in Fig. 2a illustrate that the friction 
level characterized by the width of the static curve 
hysteresis, and the LH voltage step can be quite 
different for different ETBs. 

(b) Aging. Fig. 2b indicates that the friction and LH 
voltages increase with aging of ETB. The LH 
position has been found to vary with aging as well. 
Aging may also decrease the DC motor voltage and 
torque constants due to the weakening of the motor 
permanent magnets. 

(c) Variations of external temperature and battery 
voltage. The DC motor armature resistance can 
vary in the range 1:2 as the consequence of 
temperature variations (Hashimoto et al., 2003). 
This has direct implications on the friction and LH 
voltage parameters, as illustrated by the static curves 
in Fig. 3. The battery voltage variations have similar 
influence on the process static curve parameters. 


In order to deal with process parameter variations, 
two control concepts can generally be applied: (i) robust 
control and (ii) adaptive control. The former concept 


assumes that the process model uncertainties can be 
captured by means of a specific control structure. A 
sliding mode electronic throttle control structure is often 
proposed either in its standard form (see e.g. Hashimoto 
et al., 2003; Rossi, Tilli, & Tonielli, 2000; Yokoyama, 
Shimizu, & Okamoto, 1998) or in a more advanced 
neural network-based form (see e.g. Baric, Petrovic, & 
Perić, 2002, 2004). The effectiveness of such an 
approach can be limited by high level of noise in the 
controller output signal (chattering), and/or the require- 
ments on very low sampling time and powerful 
processor. The latter concept of adaptive control 
assumes on-line tuning of controller parameters based 
on the results of on-line process identification. Such a 
concept has not been well investigated in the available 
literature, with an exception of Hashimoto et al. (2003), 
and Gagner and Bondesson (2000), where adaptive LH 
compensators are presented. 

This paper proposes a more comprehensive adaptive 
electronic throttle control strategy compared to those 
given in Hashimoto et al. (2003), and Gagner and 
Bondesson (2000), and presents the results of thorough 
experimental tests. The adaptation mechanism tunes the 
control strategy from Deur et al. (2004) with respect to 
wide-range variations of all process parameters. The 
overall adaptive control strategy can be implemented on 
a low-cost automotive microcontroller system with 
integer arithmetic only and typical electronic throttle 
sampling time in the range from 2 to 5ms. The 
adaptation strategy consists of auto-tuning and self- 
tuning algorithms (Fig. 1b). 

The auto-tuner tunes the controller parameters for 
each single ETB in the stage of vehicle production. It 
can also be executed during specific vehicle drive modes 
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Fig. 2. Variations of process static curves for different ETBs (a), and 
variations of process static curve parameters due to aging (b). 


(e.g. during maintenance intervals, or each time the 
engine is turned off), thus adapting the control strategy 
with respect to slow variations of process parameters 
(e.g. friction and LH variations). The self-tuning 
algorithm performs real-time estimation of those pro- 
cess parameters which can vary within a single engine 
run (e.g. armature resistance), and updates the relevant 
control strategy parameters. 

The paper is organized in seven sections. Section 2 
describes the process model. Section 3 presents the 
nonlinear control strategy, and gives some illustrative 
experimental results. A simulation analysis of the 
influence of process parameter variations on the 
electronic throttle performance is presented in Section 
4. The auto-tuner is described in Section 5. Section 6 


u[V] 


0) 10 20 30 40 50 60 70 80 90 
O[deg] 


Fig 3. Influence of armature resistance on process static curve 
(simulation). 


presents self-tuning algorithms which adapt the control 
strategy with respect to variations of armature resis- 
tance, battery voltage, and LH position. Different 
variants of armature resistance estimator are proposed, 
depending on the availability of armature current 
measurement and auto-tuning procedure. Both auto- 
tuning and self-tuning algorithms are verified experi- 
mentally. Concluding remarks are given in Section 7. 


2. Process model 


The process model is shown in Fig. 4 (see e.g. 
Scattolini et al., 1997). The high-bandwidth dynamics 
of the chopper is approximated by the proportional 
term with the gain Ken. The electronic throttle body is 
described by the well-known DC motor linear model 
which is extended with nonlinear models of friction and 
dual return spring. A dynamic friction model is used in 
simulations. It represents an extension of the Dahl 
friction model, and is described in detail in Deur et al. 
(2004). The load torque mz due to the throttle air mass 
flow has not been found to have significant influence on 
the electronic throttle behavior, and it is treated as an 
unknown disturbance. 

Fig. 5 shows the process static curve (cf. Fig. 2). The 
voltage parameters of the process static curve are related 
to the corresponding LH and friction torques as follows 
(cf. Figs. 4 and 5): 


K|AM .y 
2ULH = > 1 
LH Ka KiKi > ( ) 
Msc 
Usc =>. 2 
i Ken KK: ( ) 
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Fig. 4. Block diagram of electronic throttle process. u,tq,Uemf - 
commanded signal, armature voltage, and back electromotive force; i, 
- armature current; Mm, Ms, Mmf, mg - motor, return spring, friction, and 
load torque; Om, 0 - motor and throttle position; w,, - motor speed; Ken 
- chopper gain; Ka, Kı, K, - armature, torque, and voltage gain; J - 
total moment of inertia; 1/K; - gear ratio; 0zy - limp-home position; 
Mc, Ms - Coulomb and breakaway friction. 


Fig. 5. Process static curve. 


In general, the values of LH voltages Uzy+ and the 
slopes of process static curve Au/A0s above; ) and 
below,_) the LH position 074 are not equal. 

According to Pavković et al. (2003), the armature 
time constant is very small (Ta ~ 0.5ms), and it does 
not affect the linear process dynamics. Thus, the linear 
process model (Fig. 4 with my = 0 and m, = K;K10) can 
be described by the second-order lag model: 


Os) Kp 
G,(s) = = 4 , 3 
ols) u(s) 1+ p15 + aps? 63) 
with 
KenKaKı Ka J 
K == = —— = 4 
n=- KK ` Vg PTK ae 


where K, is the damping coefficient due to the back 
electromotive force Uem: 


Ka = K,K,K,. (5) 


If the return spring stiffness K, is small enough (i.e. if 
the slope of the process static curve is small), so that the 


following condition is satisfied: 


Ka>KiV/2K;J, (6) 
the linear process dynamics can be represented by a 
simple integral + lag model (IT; model) (see Deur et al., 
2004): 

Hs) K, 
uls) SC + Tems)’ 
where the process gain K, and the electromechanical 
time constant Tm are given by 


G,(s) = 


(7) 


Kp = KenKi/Ky, 
Tem =J/Ka. (8) 


According to identification results in Pavković et al. 
(2003), the condition (6) is satisfied in the region above 
the limp home position (0 > 0z, cf. slopes of the process 
static curve in Fig. 2 for 0< 0zy and 0>0zy). 


3. Nonlinear control strategy 


Fig. 6 shows the block diagram of the proposed 
electronic throttle control strategy. The control strategy 
consists of a PID feedback controller, a first-order lead- 
lag feedforward controller (FFC), and nonlinear friction 
and LH compensators. Brief description of the control 
strategy structure, tuning, and experimental verification 
is given in the next subsections. More detailed elabora- 
tion can be found in Deur et al. (2004). 


Fig. 6. Electronic throttle control strategy. 
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3.1. PID controller and feedforward controller 


Design of the linear control system with discrete-time 
PID controller is carried out in the continuous-time 
domain. The sampler and the zero-order-hold element, 
as well as the time-differentiator used in the controller 
derivative term, are approximated by a parasitic first- 
order lag term with the time constant T &Tem (T— 
sampling time). The basic set of PID controller 
parameters (for the region 0>0z7) is calculated based 
on the parameters K, and Tem of the IT; process 
model given by Eq. (7). The controller parameters 
are optimized according to the damping optimum 
analytical design method (see Naslin, 1968; Zah & 
Brandenburg, 1987; Deur, 2001). Calculating the closed- 
loop system characteristic polynomial, and equating it 
with the third-order damping optimum characteristic 
polynomial: 


A(s) = D3D3T38° + DTZ + Tes +1 (9) 


yields the following simple equations for the controller 
parameters: 


1 Tem + T 


Kr => SS; 10 

R K, DD; T? ( ) 
D:D;T. 

Tos = Dre(1 -2 a1) 

Tr =T, (12) 


By setting the characteristic ratios Də and D3 to the 
optimal value 0.5, the resulting closed-loop system is 
characterized by the quasi-aperiodic step response with 
the overshoot of 6% and the rise time of 1.87., and 
correspondingly with the gain margin, phase margin 
and critical frequency of 10dB, 60° and 2.9/T., 
respectively. The characteristic ratios Də and D3 are 
decreased here below the optimal value 0.5 (D> = 0.37 
and D;=0.4), in order to provide the boundary 
aperiodic step response of the control system. The 
equivalent time constant of the closed-loop system T 
determines the response time, and can be chosen 
arbitrarily to a value larger than 


2 T 


Pie 13 
i DD; 1 + T/Tem ( ) 


Since the return spring is significantly stiffer in the 
region 0<0@ ;y than in the region 0>6zy (Fig. 2), the 
control system tuned according to Eqs. (10)-(12) would 
have slower response in the region 0<6;7;. In order to 
provide the same control performance for both regions, 
the return spring influence must be taken into account, 
i.e. a second-order lag process model (3) must be used 
instead of the IT, model (7). The final equations for the 
optimal controller parameters for 0<0z;, can conveni- 
ently be expressed as exact modifications of the basic 


expressions (10)—(12): 


Au 
Kr = Kr — AO.’ (14) 
D3D3T? Au 
y= Tis(1 Fe Ky a) (15) 
D2D3T2 — Au\~ 
Tp_ =TH(T T)(1 a2 K, : 
D +( D+ ( TaT P a) 
(16) 


where Au/A6_ is defined in Fig. 5. The change of PID 
controller parameters with respect to operating region 
(0<0ru or 0>Oz,) is achieved by applying a simple 
gain-scheduling algorithm with bump-less transfer 
included. 

The feedforward controller introduces a discrete-time 
zero Zy into the overall closed-loop transfer function, 
thus speeding up the response with respect to the 
reference signal. The zero is determined based on the 
zero-pole canceling approach: 


zy = exp(—2T/T.). (17) 


3.2. Friction and LH compensation 


The PID controller can provide favorable electronic 
throttle behavior in the large-signal operating mode (for 
large changes of reference position Og). However, the 
electronic throttle performance significantly deteriorates 
in the small-signal operating mode due to the friction 
and LH nonlinear effects. The experimental results in 
Fig. 7a point out to unacceptable slow control system 
response with respect to reference step change of 0.2°. 
The slow response is a consequence of the stiction 
influence. Similarly, a significant response delay (a 
standstill interval) appears while the throttle passes 
through the nonlinear LH region (Fig. 8). 

In order to improve the control system performance 
in the small-signal operating mode, the PID controller 
is extended with friction and LH compensators with 
outputs uy and uzme respectively (Fig. 6). The 
compensators have the same structure, because the 
friction and LH nonlinearities have similar relay-type 
static curves and similar influences to the control 
system behavior (Fig. 4). Only one of the two 
compensators is active at any time, depending on 
whether the throttle is inside or outside the LH region. 
The LH compensator threshold parameter ezm is 
typically set to the double value of position measure- 
ment resolution Ag ~ 0.1°. 

The compensators act as bang-bang feedback con- 
trollers. Thus, they tend to provide fast and robust 
elimination of the control error ¢. The magnitudes of 
the friction and LH bang-bang controllers, Uz, and 
Usea are set near the values of the LH and friction 
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Fig. 7. Experimental step responses of systems without (a) and with 
friction compensator (b). 
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Fig. 8. Experimental ramp tracking responses of systems without and 
with friction and LH compensator. 


voltages Uzy and Us: 


Ura = KLHULH, 
Use = Ks Us, 
KLH, Ks ~ 1. (18) 


Therefore, the compensator action is limited to the 
critical nonlinear small-signal operating mode only, 
without disturbing (destabilizing) the control system in 
the large-signal operating mode. Namely, the overall 
control strategy can be regarded as a dual controller, 
where the PID controller dominantly determines the 
closed-loop system behavior in the large-signal operat- 
ing mode, and the nonlinear compensators are respon- 
sible for the small-signal operating mode. 

The main problem with the use of bang-bang 
controller is a large level of noise in the commanded 
signal u (the chattering effect), which can cause 


transmission and potentiometer wear, and excessive 
motor losses. In order to avoid the chattering effect, a 
unit-gain first-order low-pass filter is cascaded to the 
bang-bang controller (Fig. 6). In the case of friction 
compensation, the filter use is also beneficial from the 
standpoint of avoiding the step response overshoots due 
to the stiction dynamics. The filter time constants ty. and 
TrHe May be made dependent on the absolute value of 
the control error ¢ to obtain the best performance (see 
Deur et al., 2004 for explanation of tuning the lookup 
table ty-(e|)). It has been found, however, that even 
constant values ty. ~ 70 ms and try. © 15ms give good 
results for the particular electronic throttle. 

The step responses in Fig. 7 and ramp responses in 
Fig. 8 indicate that the friction and LH compensation 
results in significant reduction of stiction and LH 
standstill intervals, thus making the responses similar 
to those of the linear system. Also, the important benefit 
of friction compensation is fast establishing the one-bit 
limit-cycle steady-state response (Fig. 7). Such a 
behavior is desirable in order to decrease the control 
error even below the position measurement resolution 
(0.1° here). It is important to note that the above 
benefits are reached without any significant increase of 
the noise magnitude in the commanded signal u (cf. Fig. 
Ja and b). More detailed results of experimental 
verification of the control strategy are included in Deur 
et al. (2004), and Pavkovic (2003). 


4. Influence of process parameter variations 
4.1. Armature resistance variations 


The variations of armature resistance Ra = K7' affect 
both the parameters of linear and nonlinear parts of the 
control system (Eqs. (1)-(5)). It is reasonable to assume 
that the armature resistance variations do not exceed 
+50% of the nominal value Rao (Ra € [0.5Ra0, 1.5R,0)). 


4.1.1. Linear control system 

The linear control system consists of the PID 
controller and the linear process model given in Fig. 4, 
with my = 0 and ms = K/K,0. Fig. 9 shows simulation 
step responses of the linear control system for different 
values of armature resistance. These step responses 
point out to small sensitivity of the linear control system 
to armature resistance variations. This result has been 
proven analytically by Pavkovic and Deur (2002). The 
linear control system sensitivity is somewhat larger in 
the presence of feedforward controller (Pavkovic, 2003). 


4.1.2. Nonlinear control system 

The overall nonlinear control system is affected by 
the armature resistance change through the change of 
the breakaway voltage Us and the LH voltage Ury 
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Fig. 9. Unit step responses of linear control system without feedfor- 
ward controller for different values of armature resistance. 
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Fig. 10. Influence of armature resistance variations on performance of 
friction compensation (a) and LH compensation (b). 


(Figs. 5 and 3): 


Ka Ra 
Ury = U —=U 19 
LH LHO K, LHO Ry > ( ) 
Ka Ra 
Us = U =U: 20 
s K, 9 R (20) 


The non-adaptive control system with constant 
parameters Use = ksUso and Uppy. = KLHULHo, and 
variable process parameters Us and Uz will obviously 
have sub-optimal behavior if Ra # Rao. 

This is confirmed by the simulation step responses 
shown in Fig. 10. The step responses show that 
overcompensation of nonlinear effects occurs for R,< Rao 
(the step response is faster with possible occurrence of an 
overshoot), while undercompensation of nonlinear effects 
occurs for R,> Rao (significant standstill interval can be 
observed). The LH compensator is more sensitive to 
armature resistance variations than the friction compen- 
sator, primarily because Uz y> Us. 


4.2. Battery voltage variations 


The chopper gain K.n is proportional to the battery 
voltage Upar (Pavković et al., 2003). In the low-speed 


nonlinear operating region (Ueny * 0), the chopper gain 
Ken is cascaded to the armature gain K, = 1/R, (Fig. 4). 
Thus the influence of the battery voltage variations is the 
same in this operating mode as the influence of armature 
resistance variations (Section 4.1). 

The response in the large-signal operating mode is 
expected to be more sensitive to the battery voltage 
variations than to variations of armature resistance, 
because the gain Ken appears outside the temr damping 
path (Fig. 4). However, the battery voltage varies in a 
narrower range (typically +20% of the nominal battery 
voltage Usa, compared with +50% armature resis- 
tance variations. Thus, the influence of battery voltage 
variations to the control system response in the large- 
signal operating mode is not emphasized either, as 
illustrated by simulation results in (Pavkovic & Deur, 
2002). 


4.3. Variations of LH position 


The LH compensator is only active in a narrow region 
around the supposed LH position 0,70 + ezu (Fig. 6). If 
the information Ozypo about the LH position were 
inaccurate, the LH compensator would be activated 
outside the actual LH position 6,74. This can cause 
undesirable behavior of the electronic throttle control 
system. 

The influence of inaccurate LH position O@;y is 
illustrated by simulation step responses in Fig. 11. The 
simulation responses are compared with the first-order 
reference model response Om, which describes the desired 
behavior of the linear part of the electronic throttle 
control system. In the case when the actual limp home 
position @;y is within the LH threshold (|04 — 
OLHo| SELy), the throttle position response 0 is close to 
the model response Om, i.e. the LH standstill interval is 
relatively small and there is no overshoot in control 
system response. However, when the actual LH position 
exceeds the LH compensator threshold (0ra — 
OrH0|>éxH), a large discrepancy between the model 
response and the real response is observed, particularly 
at the actual LH position where the throttle is stuck for 
a considerable time. An undesirable overshoot can also 


lôru -0rHo |= 9.19 < Ezy — -— | Op -Orno |= 0.6° > Ery 


[deg] 


t[s] 


Fig. 11. Influence of LH position variations. 
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occur due to the inappropriate action of the LH 
compensator (Fig. 11). 


5. Automatic tuning of control strategy 


The electronic throttle control strategy in Section 3 is 
extended with an automatic tuning (auto-tuning) 
procedure, in order to deal with slow variations of 
process parameters. 


5.1. Auto-tuner outline 


The schematic in Fig. 12 outlines how the auto-tuner 
works. There are six characteristic phases of auto-tuner 
operation, which are denoted in Fig. 12 by encircled 
numbers and arrows. Each phase relates to either a 
single point or a portion of the process static curve, and 
can be executed in open loop (ol) or closed loop (cl). Fig. 
12 also shows which type of the commanded signal u or 
the reference signal Or is applied in each phase. 

The auto-tuner determines the LH position zy in the 
initial phase 0 when u = 0 is set. In phase 1 the throttle is 
positioned to the upper edge of the static curve at 0 = 
Ory, i.e. to the left most operating point of the 
approximately linear operating region 0>0,47 A 0>0 


k" phase of tuning 


open-loop/ closed loop 
operation 


shape of u (or 6p) 


0 - Limp-home position estimation 

1 - Breakaway from limp-home position 

2 - Identification of linear process dynamics 

3 - Basic tuning of controller and closing control loop 
4 - Identification of process static curve 

5 - Final tuning of control strategy 


Fig. 12. Illustration of auto-tuning operations. 


of the static curve (the nonlinear friction and LH 
influences are overcome). At the beginning of phase 2, a 
step change of the commanded signal u is applied. Based 
on the throttle position response the parameters of 
dominant linear process dynamics are estimated. Using 
the estimated process parameters, the basic set of PID 
controller parameters (for 0>677;) is determined at the 
beginning of phase 3. The throttle is then brought 
(under the PID control) to the position 69, which is 
somewhat larger than the LH position @z#. In phase 4, 
the ramp change of the throttle position reference Opr 
through the static curve region 0< 0o is commanded. In 
this way, the main parameters of the process static curve 
are estimated. Based on these estimates, final tuning of 
the control strategy is carried out in phase 5. 


5.2. Process identification 


The auto-tuning operation during the identification 
phases 0, 1, 2, and 4 is described. The description 
is illustrated by the experimental response shown in 
Fig. 13. The overall auto-tuning procedure lasts 
about 1.5s. 


5.2.1. Phase 0—LH position estimation 

The LH position zy is simply determined as the 
average throttle position during the first several instants 
of auto-tuner execution (when u = 0). 


5.2.2. Phase 1—Breakaway from LH position 

The main task of the auto-tuner during phase 1 is to 
bring the throttle to the upper point of the LH region 
(Fig. 12). This needs to be done in a relatively short 
period, using the open-loop control only, and having no 
prior knowledge of the process static curve parameters. 
Phase 1 is divided in two stages. In the first stage, the 
commanded signal u is ramped-up until the breakaway 


~ fs] 


Fig. 13. Experimental responses of electronic throttle during auto- 
tuning procedure. 
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from the LH position is detected (the interval 
t = 0-0.35s in Fig. 13). The reached commanded signal 
u is saved as the breakaway voltage upa. In order to 
avoid slow throttle transition through the region 0>0r4 
and possible excess of the upper stop position, the 
throttle is returned to the LH position in the second 
stage by setting u = 0. It is then quickly brought close to 
the breakaway operating point up, (but without break- 
away) by applying an exponential form of the com- 
manded signal u (the interval 0.35—-0.65s in Fig. 13). 

There are two practical aspects of auto-tuner im- 
plementation for phase 1 (Deur & Pavković, 2002): (i) 
choice of the rate of change of the commanded signal u 
during the first stage, and (ii) breakaway detection. The 
u-rate is chosen to a maximum value that still provides 
an accurate estimation of upa. The breakaway instant kpa 
is found as the first sampling instant after which the sum 
of three consecutive throttle speed samples is greater or 
equal to 54o (4o = 4o/T is the speed measurement 
quantization level), i.e. when œ(ks+ 1)+ w(ks +2) + 
o(ks + 3)2 54o. Using the throttle speed reconstruction 
equation w(k) = [0(k) — 0(k — 1)]/T yields the break- 
away condition 


O(ks + 3) — (ks) > 54o. (21) 
The breakaway voltage is given by 
Uba = u(Kha)s Koa = min{k;}. (22) 


5.2.3. Phase 2—Identification of linear process dynamics 

According to the analysis in Section 2, the linear 
process model may be approximated in the region 
0>0LH by (cf. Eq. (7)): 


o(s) As) K, 
uls) z uls) 1+ Toms’ 


Hence, the throttle speed step response can simply be 
described as 


Ot) = Os (1 — e™/ Fen), (24) 


The throttle speed response in Fig. 13 (phase 2) is indeed 
in agreement with the description (23). The process gain 
K, is estimated as K, = @ss/Au, where @ss is the speed 
response steady-state value, and Au is the commanded 
signal step change. In order to avoid time-consuming 
numerical optimization of Tem (applied by Pavković et 
al., 2003 for off-line process identification), it is 
convenient to integrate the speed response and find a 
relation between the speed integral (i.e. the displacement 
Ofin = Dinit) and Tem: 


T fin 
f œdt = Ofin — Dinit 
0 


T fin 
= os | d- et/Ten) dt 
0 


= Wss[T fin — Tem(1 — e77 Tomy), (25) 


Gyo() = 


(23) 


Since the identification interval Tj, is larger than 3Tem 
(the steady-state condition), the exponential term on the 
rightmost side of Eq. (25) is less then 0.05 and may be 
neglected. Hence, 


Ofin E Oinit 
Wss 


Tem x T fin _ (26) 

In order to calculate the steady-state speed w,, and 
determine Tan, it is necessary to detect the steady-state 
portion of the speed response. The steady-state instant 
kss is determined as the first sampling instant after which 
the absolute value of the sum of N consecutive throttle 
speed differences is less than or equal to the speed 
quantization level Ao, i.e. 


kss = min{k;}, (27) 
where k, satisfies the relation 


N 
XO olk, +i) — olk, + i= 1)]< Ao. (28) 


i=1 


The final throttle position Ofn and the steady-state speed 
(ss are then determined as 


Ofin = O(Kss), (29) 
1 N+M 

Oys ae @ kss +i > 30 

NFM K i i Am 


where N, < Nis the number of last several samples inside 
the sampling window used in Eq. (28), and M, is the 
number of additional subsequent samples. Typical 
values of these constants are N = 5, N; = 3, M=7. 


5.2.4. Phase 4—Identification of process static curve 

The required process static curve parameters Us, 
Urm, and Au/A@_ (Fig. 5) can be reconstructed 
from the relatively narrow portion of the overall 
static curve around the region <0 in Fig. 12. The 
static curve is recorded as the closed-loop response 
with respect to ramp form of the position reference Oz. 
The ramp slope is set to a maximum value that still 
provides a reasonably accurate prediction of the static 
curve. The recorded static curve is shown in Fig. 14 
(solid line). 

Estimation of the static curve parameters by Pavković 
et al., 2003 was based on off-line straight-line interpola- 
tions of the static curve. This technique is not suitable 
for on-line estimation, because of the high computa- 
tional requirements and certain inaccuracy of the on-line 
recorded static curve (cf. static curves in Fig. 14). 
Therefore, another identification method based on the 
characteristic points of the static curve is proposed. 
There are six characteristic points (Fig. 14), which can 
be divided in two characteristic groups: the breakaway 
points BA1-BA3, and the sliding-regime points S1-S3. 
The breakaway point BA3 was already identified in 
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Fig. 14. Process static curves recorded by off-line and on-line 
experiments, and results of on-line identification, denoted by ‘*’ and 
‘o’ symbols. 


phase 1 (ug43 = uba), and the other two breakaway 
points BAI and BA2 are identified in the similar way 
(see Deur & Pavkovic, 2002). The voltage levels of the 
sliding-regime points S, j= 1,2,3 are estimated by 
averaging the commanded signal u over a narrow sliding 
interval: 


L 


1 
ug => 5 u(k s; + i), (31) 
= 


where L is typically 5-8 samples. 

Averaging starts at the sampling instant ks; which 
is determined based on the following conditions 
(Fig. 12): 


Or(Ksj) < OLH for j = 1, 

OR({Ksj) > Or(ksj — 1) for j = 2, 

Or(ks;)>0LH A Upar already estimated for j = 3. (32) 
Once the characteristic voltage levels ug4; and us; are 
estimated (see * and o-marks in Fig. 14), the required 


static curve parameters Us, Uzm, and Au/A@_ (Fig. 5) 
are readily calculated as 


liu —u 1 /us3 — u u —u 
U= e Si (= BAI j BA2 2), 


2 2 2 2 2 
(33) 
1 /uga3 + uUsı Ugal + Us3 

= 4 
Uru A 2 2 ). (34) 

Au l fus3 —Upar uga — Hs) 
= + ; 35 
A0_ 2 (= = Omin OLH aR Omin ( 


5.3. Calculation of control strategy parameters 


5.3.1. Phase 3—Calculation of basic set of controller 
parameters 

The basic set of PID controller parameters (for the 
region 0 > 677) is calculated according to Eqs. (10)—(13). 


5.3.2. Phase 5—Calculation of other control strategy 
parameters 

The optimal PID controller parameters for the region 
below the LH position (0< 6 zy) are calculated by using 
Eqs. (14)-(16) from the basic set of controller para- 
meters (Phase 3) and the estimated slope of the process 
static curve below the LH position, Au/A0_. The 
parameter 0; of the control strategy (Fig. 6) is set to 
the LH position value estimated in phase 0. The 
parameters Use and U;;,. of the nonlinear friction and 
LH compensators are calculated based on Eq. (18) and 
estimated values of friction voltage Us and LH voltage 
Ury. 


5.4. Experimental verification of tuned control strategy 


In order to test the tuned control strategy, the auto- 
tuning procedure was extended with an interval of 
control strategy verification. The control system re- 
sponse during this interval is shown in Fig. 15. The 
control system step response has the desired aperiodic 
form for all the tested operating modes, and the settling 
time is approximately 70 ms. The estimated parameters 
were found to agree well with the results of off-line 
identification (see e.g. Fig. 14). It should be noted that 
numerous auto-tuning experiments have been executed, 
and the estimated process parameters have always taken 
on similar values. This points out to the robustness of 
the auto-tuning strategy. 


60 
40 


20 


1.8 2 2.2 24 2.6 2.8 3 
ts] 


Fig. 15. Experimental response of electronic throttle control system 
tuned by auto-tuner. 
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6. Self-tuning of control strategy 


In order to deal with the variations of battery voltage, 
armature resistance, and the LH position, which can 
occur within a single engine run, the control strategy is 
extended with a self-tuning procedure. It performs on- 
line estimation of process parameters and on-line 
calculation of the control strategy parameters. 


6.1. Adaptation to variations of battery voltage 


Since the battery voltage Upar is measured in produc- 
tion vehicles, and Ken is proportional to the battery 
voltage Ubar its variations can be compensated for by 
scaling the controller output u (cf. Fig. 4): 


us = UU pao / U bat» (36) 


where Ubaro and Upar represent nominal (initial) battery 
voltage and measured battery voltage, respectively. 


6.2. Adaptation to armature resistance variations 


The optimal parameters of the friction and LH 
compensators U;y, and Use for any value of armature 
resistance R, are calculated according to Eqs. (18)—(20), 
provided that the initial values of process static curve 
parameters Uzao and Uso are known, and that the 
relative armature resistance R4/Rao can be estimated on- 
line. The initial values Uzzo and Uso are obtained by the 
control strategy auto-tuning (Section 5). If the auto- 
tuner would not be executed, only partial adaptation of 
the nonlinear control strategy is carried out. 


6.2.1. Direct estimation of Ra/ Rao 

If the armature current measurement is available, the 
armature resistance can easily be estimated as the ratio 
of mean armature voltage and mean armature current 
for the throttle at standstill (w,, ~% 0; Fig. 4) 


K ch 


a Ug 
Ra == 


la 


(37) 


Omn~0 la 


nO 


INITIAL 
AUTO-TUNING 


Mean values are used to suppress the noise in the 
current and voltage signals i, and u. 

The block diagram of the proposed self-tuning 
strategy is shown in Fig. 16. The standstill condition 
Mm % 0 is checked by monitoring the throttle speed 
signal œ (already available from the PID controller 
derivative term), with throttle zero speed threshold Awo 
set to double value of speed measurement resolution 4o. 
Additional conditions on armature resistance estimation 
are (i) the mean armature current is not close to zero 
(to avoid division by zero in Eq. (37)), and (ii) 
the throttle is outside the highly nonlinear LH region 
(10 — 0LH| > £L). ` 

The estimated signal R,/R,o will include some 
noise and even possible spikes due to the noise in the 
estimator input signals 7, and u, and “non-ideal” enable 
logic. In order to suppress this noise, the estimated 
signal R,/Ra is filtered by a unit-gain low-pass filter 
(Fig. 16), whose bandwidth is chosen empirically to 
approximately 0.2rad/s as a trade-off between fast 
estimation and low noise sensitivity. 


6.2.2. Indirect estimation of R,/ Rao 

Simulation and experimental results in Deur et al. 
(2004) have shown that the steady-state point of the 
closed-loop control system falls on the middle curve of 
the hysteretic process static characteristic (so-called 
zero-friction curve, dotted curve in Fig. 5): 


aus 
A04 


|uss+| = ULH+ + |0 — Ory. (38) 


Since the process static curve is proportional to the 
armature resistance (Eqs. (19) and (20)), the steady-state 
controller output is rewritten as 


Au 
Ado+ 


Ra 
Rao 


| Uss+ | = 


Ra 
(Vims + 0 — Ozu) = >— [usso], 


Rao 


Fig. 16. Block diagram of self-tuning to armature resistance variations with available armature current measurement. 
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where |usso+| is the steady-state controller output 
calculated from initial process parameters U;70+ and 
Au/AOo+- 


According to Eq. (39), the ratio R,/R, can be 
indirectly estimated as 


Ra _ |Uss+| 
Rao |Usso-+ 


(40) 


The block diagram of the proposed self-tuning law is 
shown in Fig. 17a. The steady-state controller output for 
initial process parameters, |u,,9|, is calculated according 
to Eq. (39). The switch SW in Fig. 17a transfers the ratio 
|u| /|tesso+| if 0>0rH, or u| /|Usso—| if 0<@zy to the 
smoothing low-pass filter whose output represents the 
estimated ratio |uss+|/|ussot|. The filtering is enabled if 
the control error ¢ is small (the quasi-steady-state 
condition, |é|<A,), and if the throttle position 0 and 
the throttle position reference Og are outside the highly 
nonlinear LH region. The threshold parameter A, is 


(a) 


(b) 


typically set to double value of position measurement 
resolution Ag. 


6.2.3. Estimation of LH voltage Ury only 

If the initial auto-tuning procedure is not available, 
only the LH voltage Uzy can be estimated. The 
expression for the LH voltage, given by Eq. (38), can 
be simplified to 


> Au 
5 
AO, 


ULH+ ca |uss+l, if |Uss+| |0 = Ox, vO. (41) 


According to Section 2, the condition on |uss+| in Eq. 
(41) is not satisfied well in the region 0<6,7, due to the 
return spring influence. Since this may introduce 
substantial overestimation of the LH voltage for 
0<0zry, it would result in the LH overcompensation 
effect. In order to avoid the overcompensation effect, it 
is recommended that the LH correction factor kr used 
in Eq. (18) is decreased for 0< 0z. 


| Uss+ | 


low-pass filter 


Fig. 17. Block diagrams of self-tuning to armature resistance variations without armature current measurement: with initial auto-tuning (a) and 


without initial auto-tuning (b). 
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The block diagram of the proposed self-tuning 
strategy with respect to LH voltage variations only is 
shown in Fig. 17b. The LH voltage is estimated as the 
filtered absolute value of the controller output (cf. Eq. 
(41)), where a low-pass filter is applied to suppress the 
noise in the estimated signal. The enable logic is the 
same as with the adaptive scheme in Fig. 17a. The 
correction factor Kz is modified by utilizing a lag term 
with the relay input, with correction factors KLH+ 
chosen empirically (Pavkovic & Deur, 2002). 


6.3. Adaptation to LH position variations 


The analysis in Section 4.3 may suggest that the 
negative influences of LH position variations could be 
avoided by widening the LH threshold zy. However, in 
this way the friction compensator would be disabled in a 
relatively wide region around the assumed LH position 
(cf. Fig. 6), thus resulting in poor friction compensation 
performance in the region 0,74 + ery (Pavković & Deur, 
2002). Therefore, an adaptation algorithm based on on- 
line LH position estimation is proposed below. 

The block diagram of an on-line LH position 
estimator is shown in Fig. 18. The estimation is enabled 
if the throttle is stuck (|@|< A, 0) in the anticipated LH 
region (|0 — 0LHol| < AoLH), and if the model following 
error is larger than a threshold ({é,,|>A,). In order to 
avoid the overshoot effect illustrated in Fig. 11, and 
related false detection of the LH position, the enable 
logic is extended with the condition sgn(0,,, — 0) = 
sgn(Or — Om), or (OR — Om)(Om — 0) >0 (as implemented 
in Fig. 18). The threshold parameters of the LH position 
estimator are set to the following values: Agory œ% 54ọ, 
o0 = 0.14, and A,<2U,4T1/(NT KR) (see Pavković, 
2003). 

If all of the above conditions are satisfied, the enable 
signal in Fig. 18 is in the high state, and the current 
throttle position is held as the LH position 0ra 
(sample & hold block in Fig. 18). This value is 
filtered in order to obtain bump-less adaptation of the 
control strategy. 


reference 
model 


6.4. Real-time simulation of process parameter variations 


In order to provide verification of the proposed self- 
tuning control scheme, a real-time simulator of process 
parameters variations has been developed (Fig. 19). 
The simulator of armature resistance variations 
Au(AR,/R,o) is based on the motor armature model 
included in Fig. 4 (Pavkovic, 2003). Since the armature 
resistance change AR, is supposed to be relatively slow, 
it is simulated as the step response of a low-pass filter. 
The simulator of LH position variations is implemented 
by presetting the false LH position data 07, #0 to the 
control strategy (see Fig. 19, and Pavković, 2003). 


6.5. Experimental verification of self-tuning strategy 


In all of the experiments, a random low-amplitude 
throttle reference changing around the LH position is 
used as a kind of a worst-case reference shape (see 
Fig. 20a). The randomized reference is implemented in a 


Simulator of the armature resistance variations 


CHOPPER 


position variations 


Fig. 19. Real-time simulator of process parameters variations. 


low-pass 


filter 


Fig. 18. Block diagram of self-tuning strategy with respect to variations of LH position. 
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w/ adaptation —-—- w/o adaptation 


- actual R/R,, 
—— estimated R/R 


0 5 10 15 20 25 30 35 40 45 50 
(b) tis] 


Fig. 20. Control system responses (a) and estimator response (b) for 
direct estimation of R,/R,0- 


0 5 10 15 20 25 30 35 40 45 50 
t[s] 


Fig. 21. Estimator response for indirect estimation of R,/R,o- 


w/ adaptation —-— w/o adaptation 


(b) t[s] 


Fig. 22. Details of the throttle position responses for direct (a) and 
indirect (b) estimation of R,/R,0- 


form of the low-pass filtered band limited white noise 
with a constant offset (Pavković, 2003). 

The experimental results shown in Figs. 20-24 
illustrate the effectiveness of the proposed self-tuning 
schemes. The direct estimation of armature resistance 
(Fig. 20) is accurate, with relatively small ripple in the 
steady-state response. The response (estimation) delay 
of approximately 5s is caused by the narrow bandwidth 
of the applied low-pass filter, needed to reduce the ripple 


t[s] 


—--- actual Org 
estimated 0, 7, 


I 2 3 4 2 
(b1) tis] (b2) tis] 


Fig. 24. Estimation of zy (a), and detail of throttle position response 
before (b1) and after 0,7 is captured (b2). 


in the estimated relative armature resistance R,/R,o. 
This delay is not emphasized when compared to 
relatively large delays of throttle temperature transients 
that cause the armature resistance change. The indirect 
estimation of armature resistance is also accurate (Fig. 
21), but it is somewhat slower due to narrower 
bandwidth of the low-pass filter. The details of throttle 
position responses for direct (Fig. 22a) and indirect 
estimation (Fig. 22b) illustrate that the application of 
the adaptive control strategy significantly reduces the 
tracking error caused by the LH and stiction influences. 

Fig. 23 shows the experimental results for the 
estimation of LH voltage Uzy only. The response of 
estimated LH voltage Uz;y is characterized by a 
relatively large steady-state estimation error (overesti- 
mation). The error is compensated for in the estimator 
output signal Uzme by applying a variable compensation 
correction factor Kz4(0). 

The experimental responses for the estimation of LH 
position variations are shown in Fig. 24. The estimator 
of LH position variations correctly captures the actual 
value of LH position (Fig. 24a). The delay of the 
estimator response is caused by a conservative choice of 
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the model following error threshold A,, which is needed 
to provide robust estimation. The comparative response 
details in Fig. 24b clearly indicate the benefits of the 
application of LH adaptation: the adaptation results in 
significant reduction of the LH standstill intervals and 
throttle position overshoots. 


7. Conclusion 


A nonlinear strategy of electronic throttle control has 
been presented. The strategy includes linear feedback 
and feedforward controllers, and nonlinear friction and 
LH compensators. The presented experimental results 
have shown efficient rejection of the friction and LH 
effects that resulted in fast and accurate control system 
responses for the whole electronic throttle operating 
region. 

A simulation analysis has been carried out in order to 
investigate the influence of process parameter variations 
on the behavior of the electronic throttle control system. 
The battery voltage and armature resistance variations 
result in proportional variations of the different process 
static curve voltage parameters, which, in turn, sig- 
nificantly affect the behavior of the nonlinear part of 
control system. On the other hand, the linear part of the 
electronic throttle control system is robust to the 
armature resistance and battery voltage variations. 
The analysis of the robustness with respect to LH 
position variations has pointed out to significant 
deterioration of the control system performance if the 
LH position varies outside the LH compensator 
activation threshold. 

In order to improve the control system robustness, the 
control strategy has been extended with an adaptation 
mechanism consisting of auto-tuning and self-tuning 
algorithms. The auto-tuner application provides that all 
electronic throttles from a vehicle production series have 
similar (desired) control performance, regardless of 
variations of electronic throttle body parameters due 
to production deviations, external conditions variations, 
and aging. The auto-tuner does not require any prior 
knowledge of the process parameters. It is characterized 
by simple implementation and short execution time 
(about 1.53). 

The self-tuning strategy has been developed in order 
to deal with process parameter variations which occur 
during a single engine run (i.e. those that cannot be 
captured by occasional utilization of the auto-tuner). 
These parameters are the battery voltage, the armature 
resistance, and the LH position. Since the battery 
voltage is measured in production vehicles, the self- 
tuning adaptation to battery voltage variations can be 
readily done by means of a gain-scheduling algorithm. 
The other two parameters, the armature resistance and 
the LH position, are estimated on-line based only on the 


standard electronic throttle measurement signals. The 
self-tuning strategy has been experimentally verified for 
the critical case of slow and stochastic changes of the 
reference throttle position through the LH region. The 
experimental results have demonstrated fast and accu- 
rate process parameter estimation, and consequently 
superior performance of the adaptive control strategy 
compared to the classical, non-adaptive strategy. 


Acknowledgement 
This work has been supported by Ford Motor 


Company and the Ministry of Science and Technology 
of the Republic of Croatia. 


Appendix. On the block diagram formulation 


Single-box blocks show step responses of linear terms. 
Two examples are given below. 


K 
Proportional *¥ & ¥ 
term Y=Kpx 
First-order = y T; a y=x 
lag term dt 


Double-box blocks correspond to input-output maps of 
nonlinear terms. These are some examples: 


Sien org y =U sgn(x) 
(times 
constant) 

0, | x I< Xth 
Switching x Yo y= 
logic A 1, otherwise 
Look-up table k » y=f(x) 


+,—, X, +. 
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